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For the Grimus-Stockinger formula the same dimensionless parameter of asymptotic 
expansion is found by several ways of calculations. This parameter strongly depends on 
the width of wave packet. 

PACS: 14.60.Pq 

u?\ 1 Introduction 

For the modern theory of neutrino oscillations [H |2] the main tool is Grimus-Stockinger 
theorem [3], which gives the leading asymptotic behaviour with |R| = R — > oo for the integral: 
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where: R = Rn, n 2 = 1, and the function $(q) G C 3 decreases at least like 1/q 2 together 
with its first and second derivatives. In order to understand the physical conditions necessary 
for this expansion the dimensionless parameters should be determined from the higher order 
corrections to this formula. Here this parameter is defined unambiguously by the use of various 
asymptotic expansions allowing to calculate the further corrections. 



2 Corrections for three-dimensional case. 

To obtain the higher corrections of order R~ n we suppose that $(q) and its first and second 
derivatives are represented by Fourier-transform as: 

$(q) = fd 3 xe i(q ' x) y?(x), V g $(q) = i fd 3 x e i(qx) xy>(x), and so on. (2) 

Since 1/q 2 is also Fourier-image of l/|x|, Eqs. (j5]) are valid at least in the sense of distributions 
also for the functions defined above in [5J. By using the first equality of the following well 



known representations for spherical wave as a free Schrodinger 3- dimensional Green function 
with k = 2A, q = 2p: 
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and by interchanging the order of integration for integral ([T]) one finds: 
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Substituting here the expansion, which in the exponential should always contain one additional 
order with respect to the ones in denominator: 
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we come to the corresponding expansion of integral PJ up to 0(R 
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that by making use of (J2J) transcribes as: 
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with: (n ■ V g )$(q)| q= _ Kn = — d K Q(— kxi), and so on. 
For any positive definite quadratic form of momentum q: ( = (qA _1 q) > 0, 

with: $(q) = U((), a(n) = (nA _1 n) , c?(n) = (nA~ 2 n) , that is: 
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Then for Gaussian wave packet: %{£) = e ^ 4 , the expression (jSJ) reads: 
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Here the square bracket evidently represents corrections only to the phase of the exponential. 
It may be directly obtained by the saddle-point method. 

To this end let's transcribe the integral ([1]) for above Gaussian wave packet by using the 
second representation of Eq. (J3J). Gaussian integration gives: 



J( R ) = i 
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where: iF{t) = it(\ 2 + iO) - (RK t R) , iF"{t) = 2 ( R {K t } 3 R 



iF'{t) = i [A 2 + iO + (R {KJ 2 r)] 1-4 0, t = R/X + ia(n) + e 
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- is the saddle point as asymptotical solution of Eq. ( ITT]) for R — » 00 up to the e = 0(X/R). 
It is obtained by diagonalization A = O t AO onto the eigenvalues A = diagja.,}, with 
< a, = I/a, < 00 and determinant |A| = det{A} = aia 2 a 3 , by using a suitable orthogonal 
rotation: g = OR, g 2 = R 2 , and due to Eq. (TTTT) defines F(t ) and |K t J up to 0(e 2 ). Along 
the path deformed according to a(n) > we obtain: 
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that for: k = 2A, TrjA" 1 } = f^ 
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exactly coincides with Eq. (jSJ) with the same precision. The corrections in (jSJ), f fT3|) evidently 
disappear for degenerate case: atj = «i, for j = 2, 3. 
For the neutrino oscillations problem: k = 



E 2 . — m 2 w I£ re , and for the Gaussian wave 



packet with coordinate width o x : A ~ cr~ 2 , so Tr{A -1 } ~ a(n) ~ <r 2 , whence, the true 
expansion parameters appear as combinations of two different dimensionless ones: k<j x and 
<J X /R, that define the application conditions of Grimus-Stockinger formula as: 
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3 The four- dimensional case. 

In fact the above integral ([1]) is only three-dimensional part of the four- dimensional one defin- 
ing macroscopic Feynman diagram [2] of the problem: 
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for a 2 = iO — x 2 = e l7T x 2 , is the causal propagator in coordinate space, and now the four- 
dimensional Fourier representation is assumed for $(g), which for relativistic Gaussian wave 
packet |2] reads as: 
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Here again A~ l ~ a 2 in terms of Gaussian coordinate width for any positively defined 
quadratic form of momentum q in Minkowski space: ( = (qA~ l q) > 0. Then for a 2 — > 
one has: A — > oo, $(g) i-> 1, 0(r) i— >■ 5 4 (r), whence: iJ(R) h> m 2 h(i0 — m 2 R 2 ), that is 
reasonable from physical viewpoint. 

Repeating now all the previous steps (J!])-® for the second expression (JIB]) of J(-R), with 
the so approximated propagator (TTB1 . one obtains for arbitrary $(g) (fT9l) and 77^ = R^/vR 2 , 
l = i\l\, \l\ = y/R 2 , with 
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that for the relativistic Gaussian wave packet from (TT9I) . with 0(77) = fr;i 77), a 2 (7/) 
(t]A~ 2 7]j, analogously gives: 
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Such, on the first sight, rough calculations are exactly confirmed again by saddle point method. 
Indeed, by exponentiating like ([3]) the denominator of the first expression (1161) with m = 2m 
and representation dependent g \-¥ g^u or 5^, by means of Gaussian integration clarified in 
Appendix one has instead of (j9"l)- (fTT|) : 
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to is again the saddle point as asymptotical solution of Eq. ( 125]) for |/| — > oo up to the 
e = 0(m/\l\). It is obtained now by diagonalization in Minkowski space as A = ^'(A)j n ^ n ' 
onto the eigenvalues (A)j n = gj n <i{j) ordered jlj as < gjjd(j) = gjj/oij < oo, with determinant 
\A\ = \A^ U \ = \(A)j n \, by using a suitable Lorentz transformation as gi = ^^-R^ with 
8j = 9jj^ = 9jkQ k , Q 2 = R 2 , and due to Eq. (|25|) defines again J-"(to) and |JC(t )| up to 0(e 2 ). 
Along the path deformed according to airf) > 0, instead of (TT2|) one finds: 
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exactly coincides with Eq. fT22|) with the same precision. 

The true parameters of expansion appear again as the following products of the two di- 
mensionless parameters, that are now ma x and a x /\l\: 



{ma x )-rjj < 1, and: {ma x f -y < 1, 
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and they have the same order for (ma x ) < 1. 

It is easy to see that both the conditions for 3- and 4-dimensional cases are practically the 
same. Indeed, Eq. ( TT71) implies that 
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Since for ultrarelativistic neutrino T 
rewritten as: 
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Thus, for (ma x ) < 1 these both conditions are the same as the first one in 3 - dimensional 
case (fl5]) . Moreover the same dimensionless parameter (1321) defines in fact the asymptotical 
solutions of both the saddle points equations ( ITT]) and (T25J) . Note that exact values of the 
first's and second's square brackets in (fT3|) and/or (T271) respectively may be different, and 
their determination in terms of a x for 4 - dimensional case [2] ([2"T1) - ([2"9"1) is different from 3 - 
dimensional case (fTB"j) . flHJ). 

The authors thank V. Naumov, D. Naumov, E. Akhmedov, S. Lovtsov, and N. Iljin for 
useful discussions. 



4 Appendix. 

In order to strictly calculate a standard Gaussian integral over the Minkowski space [2J: 

Jd%exp{-(yAy) + 2(By)}, (34) 

where the quadratic form (yAy) = y^A^ u y u is symmetric and positive definite, the following 
solution of the eigenvalue problem may be used: 
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In spite of ambiguity ([5] §94) of diagonalization procedure for symmetric tensor in Minkowski 
space, the positive definiteness of A leaves the used type of its diagonalization only, leading 
to eigenvectors ^ n \ n = -j- 3 (1351) . whose components ^ n - ) define Lorentz transformation 
diagonalizing the form. Then with the substitutions Y n = £,^y p transforming {yAy) = 
(YAY), b m = B^\ (By) = (bY) = b m g mn Y n , the integration (|33) factorizes to: 
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where A^ (A" 1 ) = <J/, and \A 



^ , — ^ i^x^i = I (A) fa | is denned by 
Nevertheless it is instructive obtain the same result without reference to diagonalization by 
using the direct integration over space and time variables separately. Since for n-dimensional 
Minkowski space with signature metric g^ = diag{l, — 1, — 1, . . . , —1} in any given orthogonal 
basis the symmetric tensor A is represented by the block matrix A M „ = Aij for i,j,fi,v = 
-r n — 1, with the rightmost bottom block A%j = A 1 -* = Aji, for i, j — 1 -r- n — 1: 
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the integral (1341) with d 4 y i— >■ Gpy may be rewritten as: 
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The both integrals are over Euclidian space now, so they are evaluated to 
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if a = A 00 — A ; (A -1 ) A fe0 > 0. The expression is simplified by Laplace expansion of the 
determinant |A^|, where M Q* ^ '") m eans the minor of the matrix A^, whose rows ii, ii-, ■■■ 
and columns ji,J2, ■■■ are deleted: 
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Thus a > due to positivity condition of the form ( 1371) implying that |A^|, \Aij\ > 0. 
Furthermore, if for any symmetric block matrix A: 
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whence the rightmost bottom block of the inverse to (j37p is expressed for i, k, l,j = 1 -r n— 1 



as: 
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the argument of the exponential in (J55|) is also reduced to expression (155|) : 
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A generalization of the integral (13"4"j) with arbitrary polynomial or smooth function similar 
to the well-known approximations for Euclidian case [B] here is also straightforward. The 
imaginary part of A (1371) in integral ( 134")) appearing for (1231 as itg^ preserves convergence of 
the integral and may be included by analytical continuation of the expression ( 136]) . 
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